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Abstract 

In a multi-user wireless network equipped with multiple relay nodes, some relays are more 
intelligent than other relay nodes. The intelligent relays are able to gather channel state information, 
perform linear processing and forward signals whereas the dumb relays is only able to serve as 
amplifiers. As the dumb relays are oblivious to the source and destination nodes, the wireless network 
can be modeled as a relay network with smart instantaneous relay only: the signals of source- 
destination arrive at the same time as source-relay-destination. Recently, instantaneous relaying is 
shown to improve the degrees-of-freedom of the network as compared to classical cut-set bound. 
In this paper, we study an achievable rate region and its boundary of the instantaneous interference 
relay channel in the scenario of (a) uninformed non-cooperative source-destination nodes (source 
and destination nodes are not aware of the existence of the relay and are non-cooperative) and 
(b) informed and cooperative source-destination nodes. Further, we examine the performance of 
interference neutralization: a relay strategy which is able to cancel interference signals at each 
destination node in the air We observe that interference neutraUzation, although promise to achieve 
desired degrees-of-freedom, may not be feasible if relay has limited power Simulation results show 
that the optimal relay strategies improve the achievable rate region and provide better user-fairness 
in both uninformed non-cooperative and informed cooperative scenario^. 

Index Terms 

This work has been performed in the framework of the European research project SAPHYRE, which is partly funded by 
the European Union under its FP7 ICT Objective 1.1 - The Network of the Future. 
'Part of this work is submitted to ISIT 2012 [1]. 
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I. Introduction 

In a wireless network, where the source (S) nodes and destination (D) nodes communicate in a 
shared medium, e.g. time, frequency, code, the signals interfere with each other in the air and may 
not be possible to be separated at the destination nodes which then leads to a rate degradation. 
Interference management techniques are thereby crucial to the ever-increasing demand of quality 
of service. Advanced interference management techniques, e.g., superposition coding and multi-user 
decoding, induce difficulty in practical implementation as all codebooks and interference decoding 
capabilities are required respectively. We assume in this work that the S and D nodes are not able 
to perform multi-user encoding and decoding techniquesP. To improve the rate performance of the 
system, one can introduce relays to the system and therefore obtain an interference relay channel 
(IRC). The relay is responsible for signal boosting and interference managing. We assume that the 
relay employs an amplify-and-forward (AF) strategy which provides flexibility in implementation as 
the relay is oblivious to the modulation and coding schemes in the communication between S and D 
nodes [3]. 

The novel notion of relay-without-delays, also known as instantaneous relays, originally proposed 
in [4], refers to a type of relays that can forward signals consisting of both the current symbol and 
the symbols in the past, instead of only the past symbols as in conventional relayj^. The authors 
further point out that a relay-without-delays is only a conventional relay with one extra delay on the 
link from S to D. A visual example is given later in [5] where in a network with multiple S-D pairs 
and conventional relay nodes, some naive relays only AF signals and some smart relays can gather 
channel state information (CSI) and forward linear processing of signals. As the naive relays do not 
affect the optimization of the system and are oblivious to S and D, the links from S to D through 
the naive relays are equivalent to direct links from S to D with one extra delay. Hence, the network 
is equivalent to an IRC equipped with smart relays-without-delays only. It is shown in [6], [7] that 
instantaneous IRC can achieve degrees-of-freedom (DOF), or the multiplexing-gain, higher than the 
classical cut-set bound but can greatly simplify the interference alignment scheme [8]. 

Interference neutralization (IN) is a technique of canceling, zero-forcing or neutralizing interference 
signals by a careful selection of forwarding strategies when the signal travel through relay nodes before 

^For the results on instantaneous relay capacity with multi-user decoding destination nodes, please refer to [2]. 
^In the terminology of half-duplex and full-duplex, instantaneous relay is a kind of full-duplex relays. 
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reaching the destination. This general idea has been applied in deterministic channels [9], [10] and in 
two-hop relay channels [3], [11], [12], also known as multi-user zero-forcing and orthogonalize-and- 
forward. To distinguish from IN, aligned interference neutralization (AIN), which is a combination 
of interference aligimient and interference neutralization techniques [5], [13], is first proposed in a 
conventional full-duplex two-S, two-relay and two-D network (the 2x2x2 chaimel) [13] and is later 
extended to the IRC with an instantaneous relay [5]. Agreeing with [6], the achievable DOF with 
an instantaneous relay is higher than the cut-set bound, achieved with the utilization of AIN. For 
example, the DOF of a two-user IRC with one instantaneous relay is shown to be 1 .5 as compared 
to 1 in a two-user interference channel (IC) without any relay. 

Although instantaneous relaying can increase the DOF of the underlying interference channel, the 
DOF analysis describes only the high SNR performance of the system, leaving the system behavior in 
low to medium SNR regime uimientioned. For any SNR regimes, the boundary of an achievable rate 
region, the so-called Pareto boundary, holds particular importance because it is highly desirable in 
the system's interest to operate on the Pareto boundary, i.e., Pareto efficient: no user can improve its 
own rate without decreasing other users rate [14]-[16]. It is the purpose of this paper to characterize 
the Pareto boundary of the instantaneous IRC with (a) the optimal AF relay processing matrix and 
(b) interference neutralization. Our work differs from the literature in the following areas. We employ 
IN in an instantaneous relay channel with direct links between S and D, whereas [3], [11], [12] 
consider relays-with-delays and with no direct links between S and D. To maintain simple design 
of S and D, we consider no symbol extensions and utilize a memory-less instantaneous relay and 
IN, whereas AIN (a combination of interference aUgnment and IN by using instantaneous relay with 
memory and symbol extensions) is employed in [5], [13]. Note that symbol extensions are shown to 
achieve desirable DOF. For details, please refer to [17]. We assume Unear relay processing which is 
considered in all aforementioned work. For non-linear instantaneous relaying, please refer to [18]. 

However, the performance analysis of the IRC is not Umited to IN. Other interference management 
techniques, such as interference ahgnment and interference forwarding, are shown to achieve promi- 
sing results. In [19], on a quasi-static channel, interference alignment is implemented on a M-user 
IRC with each node equipped with a single antenna. The power scaling of the relay is computed 
in order to guarantee a degrees-of-freedom ^ in the channel. In [20]-[22], instead of neutrahzing 
interference signals at the destination nodes, the relay is responsible for forwarding interference signals 
to the destination nodes so that the strength of the interference signal is increased and is able to be 
decoded and subtracted. However, a unified algorithm which is capable of adapting different relay 
strategies, such as interference aligimient, IN and interference forwarding, depending on the chaimel 
qualities, is yet an open problem. For the reference of the readers, we summarize the achievable DOF 
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TABLE I 

The summary of achievable DOF in different topologies of the interference relay channel. 
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Method 
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3M/4 


IR 


AIN [5], potent relay 


SUD S-D 


2 
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2 


2 


CR 


AIN [13], potent relay 




K 


1 


K 


1 


K 


IR 


IN, non-potent relay (this paper) 




K 


1 


K 


1 


K 


CR 


IN, ME, non-potent relay (Pr ~ 0(P^)) [23] 


MUD S-D 


2 


1 


(1,2) 


1 


2 


IR 


MUD, non-potent relay [2] 



of IRC with instantaneous relays and relay s-with-delays (conventional relays) in Table |£. Note that 
the degrees-of-freedom is only a performance measure in the high SNR regime and is not suitable to 
be used as a performance metric in the low and medium SNR regime. 

In this paper we study the performance, in terms of achievable rate regions and rate fairness, 
of a K users instantaneous IRC. We assume a power constraint at the relay and linear processing 
capabilities on the relay and S-D pairs. Global CSI is assumed at the relay which is a common 
assumption in aforementioned works. Although the effort of CSI estimation and processing at the 
relays may seem large, the resulted gain largely outweighs the cost in low mobility environments [3]. 
A distributed smart relays system, which is capable of gathering CSI, matching the LTE standard, is 
recently demonstrated [24]. In Section HU we describe the channel model. An achievable rate region 
and the Pareto boundary formulation for (a) AF strategy optimization and (b) AF and IN strategy 
are given in Section Uni In the same section, we address the necessary and sufficient condition of 
feasibility of IN in Theorem [T] in terms of channel conditions and the physical attributes of the relay 
(e.g. the number of antennas and the power at the relay). 

An improvement of rate performance is expected, as turning off the relay gives the same rate region 
as the single-input-single-output (SISO) IC. The interesting idea is to improve the rate performance by 
optimizing the relay strategy, even if the S-D nodes are uninformed of the relay 's presence and have 
their strategies unaltered (See Section\IVi. The computation of the optimal relay processing matrix, 

''The notations # S-D stand for number of source and destination nodes in the system. Ns.d is the number of antennas 
at source and destination nodes. Nr is the number of antenna at the relay. # R is the number of relay nodes in the 
system. DOF is the degrees-of-freedom achievable. All relays are full-duplex. IR represents instantaneous relays and 
CR represents conventional relays. Methods used to achieve the DOF can be aligned interference neutralization (AIN), 
interference alignment (lA). ME stands for Markov encoding and MUD stands for multi-user decoding whereas SUD 
stands for single-user decoding. A potent relay is a relay with unlimited power whereas a non-potent relay subjects to 
transmit power constraint. The notation (1,2) refers to one receive antenna and 2 transmit antennas at the relay. 
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optimal in the sense of Pareto optimality (see Section |III] for details), is not trivial as it closely 
resembles an IC which is well known for the NP-hard complexity of the transmit optimization [25], 
in terms of the number of users in the system. However, we are able to compute an upper bound 
using the Semi-Definite-Relaxation (SDR) techniques which relaxes the optimization problem to a 
convex problem and can be solved efficiently. As we will see later in Section |IVl in the scenario with 
two S and two D nodes, this relaxation is shown to be tight. In Section |Vl we investigate the scenario 
when the S-D nodes are aware of the relay's presence and are willing to cooperate by optimizing 
their transmit power. The resulting rate region can be computed using SDR techniques and can be 
compared fairly to the achievable rate region of a SISO IC, without any relay, in which the S nodes 
optimize their transmit power. We show numerical evidence and discussion of the advantages of relay 
in a system of uninformed S-D nodes in Section |Vll Conclusions and future directions are given in 
Section Ivni 

A. Notations 

The set C"^* denotes a set of complex matrices of size a by & and is shortened to when 
a = b. The notation M{A) is the null space of A. The set R denotes all real numbers whereas 
M+ denotes the set of positive numbers. The operator represents expectation over the random 
variable a. The operator denotes the Kronecker product. The superscripts ^ , represent 
transpose, Hermitian transpose, inverse and transpose, inverse and Hermitian transpose and Moore- 
Penrose inverse respectively whereas the superscript * denotes the conjugation operation. The scalar 
and vector Euclidean norms are written as |.| and ||.||. The trace of matrix A is denoted as tr(A). 
Vectorization stacks the columns of a matrix A to form a long column vector denoted as vec(A). 
The identity and zero matrices of dimension K are written as Ik and Oa'- The vector ej represents 
a column vector with zero elements everywhere and one at the i-th position. The notation [A]mi 
denotes the m-th row and l-th column element of the matrix A. The notation p^.^, < a < b < n, 
denotes a vector which has elements \pa,Pa+i, ■ ■ ■ ,Pb] where p = [pi, . . . ,p„]. 

II. Channel Model 

In a SISO IRC, an example of two sources and two destinations shown in Fig. [TJ we denote the 
sources as Si and destinations as Dj, i = 1, . . . , K. The multi-antenna relay node is denoted as R. 
Denote the complex channel from Si to Dj as hji and the complex channel vector from Si to R as g^j 
and from R to Dj as g^y. All channels are assumed to be independent identically distributed complex 
Gaussian variates, g^j,gjY G C^^^^, where M is the number of antennas at the relay. We assume a 
memory less instantaneous relay [2], [4], [5] which has a linear processing matrix R € C*^^^^. The 
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(a) a relay-assisted network (b) instantaneous relay network 

Fig. 1. The wireless relay-assisted network with layer one repeaters and one smart relay is shown in subfigure (a). The dotted 
hnes demonstrate the equivalent links between source the destination taking into account of the presence of the repeaters. 

All paths from source to destination nodes take two time slots and links from source to relay and relay to destination take 
one time slot. The equivalent channel is established in subbfigure (b) by replacing the relay as an instantaneous relay and 
information through instantaneous relay arrive at destinations the same time as the direct links. 



signal received at is: 

K 

yr = ^ErjX3+nr (1) 

i=i 

where Xj are the transmit symbols from Si which is assumed to be zero mean proper Gaussian variable 
and has power constraint PJ""^, E [Ixjp] = Pj < P^°'^,j = I,. . . ,K. The noise at the relay is 
denoted as which is assumed to be independent identically distributed proper Gaussian variables 
with zero mean and unit variance. The assumption of circularity for the transmit symbols simply the 
derivation as the achievable rate is the Shannon rate. For the degrees-of-freedom and achievable rates 
improvement with improper Gaussian signaling, please refer to [26], [27] respectively. The received 
signal at Dj,j = I, . . . , K, is 

K 

Vj = ^ {hji + gj^ R Sri) xi + gfj.Knr + Uj (2) 

1=1 

For brevity, denote p = [Pi, . . . ,Pk]'^ G M;^^^. The Signal-to-Interference-and-Noise ratio at 
destination j is 

SINK, (R, p) = \hjj+gfr^grj\'Pj 

j:f=i,i^j \hji + gf. Rgw m + II gf. R p + 1 

where || gj^ R p is the amphfied noise from R to Dj. The power constraint at the relay is 

EyJtr(Ry,y,^R^)] <P--. (4) 

Note that ^n^,Xj,j=i,...,K [y^y^] = X^jLi Srj g^i power constraint is therefore rewritten 

as the following: 

tr(RQR^) <P;"'^^ (5) 
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where 

K 

Q = J]g,,g,^.P,+I. (6) 



A. Interference neutralization 

In order to neutralize interference, the following K{K — 1) equations have to be satisfied at the 
same time: 

^ + g,>R-gri = 0, i,j = l,...,K,ii^ j. (7) 

Let Gdr = [gir> g2r7 • • • ; ^Kr\ ^nd Grt = [gri, gr2> • • • > ^tk]- Denote Si = Rgrf. We have the 
interference neutralization requirement, 



G% R Grt — S 



(8) 



with 



-h2i 



S2 



-h 



K{K~l) 



SK 



(9) 



Note that S is a matrix with off-diagonal elements as the channel coefficients of the interference 
channel and diagonal elements as the optimization variables s. As we sill show later, the optimization 
can be facilitated if the optimization variable is S instead of s. This is due to the linear relationship 
between R and S. However, we must stress that only the diagonal elements of S, s, are free to be 
optimized as S is constrained to the form in Eqt. (|9]l can be rewritten to a more comprehensive 
form. We introduce a row selection matrix T which select the off-diagonal elements of S from 
the vector vec(S). For example when = 2, we have T = [0, 1, 0, 0; 0, 0, 1, 0] and Tvec(S) = 
[-/i2i, -hi2]^. We have 

Tvec(S) = -Tvec(H). (10) 

If IN is feasible (see later in Theorem [T]for feasibility issue), we can choose a relay matrix R, which 
is a function of complex coefficients s, that satisfies the IN constraint dD and achieves the following 
SINR, 

(11) 



SINRf(S,PjO 



5rj 



l|gj^.R(S)P + i 

In the following, we use the matrices R and R(S) interchangeably when we wish to emphasize 
the optimization parameter S. Please see the next lemma for the formulation of R(S). In order to 
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make sure the requirements in ([8) are not over-determined, we find the minimum number of antennas 
required in the relay such that dSJ is feasible. 

Lemma 1: A sufficient condition of IN in ([8]) on the minimum number of antennas at the relay 
is M > K. For some target signal coefficients si, . . . ,sk £ C, the relay processing matrix R that 
satisfies the interference neutralization requirement (|8} is determined by 

vec(R) = [Gj^t^G^y vec{S). (12) 
Proof: From dS), using vec(A B C) = (C^ (g) A) vec(A) [28], we obtain {Gjt «> G^) vec(R) = 
vec(S). For a consistent system [29, P.12, Section 0.4.2], a sufficient condition is to have the number 
of equations (length of vec(S)) less than or equal to the number of variables (length of vec(R)) 
M > K. Multiply both sides with the Moore-Penrose inverse of (G^^ (g) G^) and result follows. ■ 
Note that when M < K, the system in ([8]l has more equations than unknowns and the consistency 
of the system depends heavily on the particular channel realizations. For simplicity of presentation, 
from now on, we set M = K. We assume that the matrix (G^® G^) is full rank and obtain [28, 
P. 35] 

vec(R) = (G^t®G^,)"%ec(S). (13) 

If there is no power constraint at the relay, given any target signal coefficients Si, we can construct 
a relay matrix R as in (IT3] ) such that the IN requirement is satisfied. Otherwise, any performance 
metrics, e.g. sum rate, subject to the IN and power constraints, are optimized in a dimension of K 
by optimizing the complex coefficients s = [si, . . . , skY . 

III. The Pareto Boundary optimization problem 

The achievable rate region of an instantaneous AF relay is defined to be the set of rate tuple 
achieved by all possible relay processing matrix R satisfying the power constraint ©: 

n= U (C(SINRi(R)),...,C7(SINRx(R))) (14) 

R:tr(RQR")<P, 

where C{x) = log2(l + x). Similarly, we define the achievable rate region of an instantaneous AF 
relay with IN to be the set of rate tuple achieved by all possible relay processing matrix R satisfying 
the IN constraint (IT3]) and the power constraint (|5): 

7^™= U (C(SINRi(R)),...,C7(SINR/^(R))). (15) 

R:tr(RQR")<P,, 
vec(R)=(G^,®G«)-'vec(S), 
Tvec(S)=-Tvec(H) 

Note that the IN requirement gives a smaller feasible set and thus T?J^ C TZ. However, the motivation 
of the study of 7^™ is two-fold: (a) intuitively when both transmit power at S and R is very large. 
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the optimal relay strategy is to neutralize interference, so as to transmit at the maximum DOF; but 
the performance of other SNR regimes is yet to be studied, (b) the characterization of TZ^^ in the 
instantaneous IRC with direct link between S-D pairs is still an open problem. The outer boundary 
of TZ (TZ^^) - the Pareto boundary (PB) of TZ (7^™) - is a set of operating points at which one user 
cannot increase its own rate without simultaneously decreasing other users rate. 

Definition 1 ( [16], [30]): A rate-tuple (ri, . . . , ri^-) is Pareto optimal if there is no other rate tuple 
{qi, ...,qK) such that {qi, ...,qK)> (n, ■ ■ ■ .tk) and {qi, ...,qK)^ (n, . . . ,ri^j£l 
By definition, the operation points on the PB can be evaluated by maximizing one user's rate while 
keeping other users' rates constant. Other optimization techniques for PB evaluation has been proposed 
[15], [31]. Here, the PB problem is formulated as the maximization of SINRi subject to the constraints 
on SINRj > 7j for some pre-determined target SINR values 7^ , j = 2, . . . , K. 

Problem Statement 1: The Pareto boundary of TZ (IT4l) is a set of rate tuple ^C(7*), C(72), . . . ,C{'~^k) 
where is the optimal objective value and 7j,j = 2,...K are the constraints of the following 
optimization problem. 



max SINRi(R,p) 



(PB)<^ 



s.t. SINRj(R,p)>7j, j = 2,...,K, (1^) 



tr(RQR^) < P,'"'^^ 

Similarly, we formulate the PB optimization problem with IN in the following. To this end, we 
combine the first two constraints in the feasibility set of (IT5]) . Using the fact that tr(ABCD) = 

vec(D^)^ (C^(g) A) vec(B) [28], we have 

tr(R Q R^) = vec(R)^ (Q"^ (g) l) vec(R) 

vec(S)^ ( (diag(p) + G,-; G,",^) G.-^ G,"^) vec(S) d^) 



Q 



-T 



rt 



where (a) is due to ©, ^ and G^* G'^^ = diag(p) + G^* G 

Problem Statement 2: The Pareto boundary of TZ^^ ([T5] ) is a set of rate tuple ^C(7*), (^(72), . . . , C(7i^) 
where 7* is the optimal objective value and 'jj,j = 2,...K are the constraints of the following opti- 

^The inequality is component-wise. 
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max SINR^iN(S,Pi) (18a) 

(PB-IN) I s-t- SINR™(S,P,) >7„ j = 2,...,K, (18b) 

vec(S)^Qvec(S) < P^"'' , (18c) 
Tvec(S) = -Tvec(H). (18d) 

mization problem. 

Before we show the optimization methods of the aforementioned problems, the feasibility issue of 
(fTSl) needs to be addressed. 

Theorem 1: A necessary and sufficient condition for the feasibility of interference neutralization - 
satisfying (I18cl ) and (I18dl ) simultaneously - is 

(Tvec(H))^ (tQ'^T^Y^ Tvec{U) < P^^"^, (19) 

where Q = f diag(p) + 0,7/ G^^) <^ (g. ^ G.^) . A feasible solution is 



vec(S) = F (x„ +(T F)t T vec(H)) (20) 

where x„ G AA(TF) and Q"^ = FF^. 

Proof: See Appendix IVIII-Al ■ 
Note that the left hand side of the condition ([T9l ) is the power of the relay matrix which only 
performs IN (with x„ = in (l20b ) and only depends on channel coefficients whereas the right 
hand side is the relay power constraint. The importance of Theorem [T] lies in the practicality of IN 
feasibility verification. With the information of channel qualities H, Grt, G^r, transmit power at S p 
and transmit power constraint at relay P^"^, we can immediately check whether IN is feasible or not. 
Further, if IN is feasible, there is always one feasible solution in (l20l ) by setting x„ = Orxi, if there 
is excess power at relay, we can optimize x„, and in turn vec(S), to improve system performance, 
which is the objective of the succeeding section. 

IV. The optimal relay strategies with uninformed S-D nodes 

In this section, we analyze the PB problems taking into consideration that the S and D nodes are 
uninformed of the presence of the relay nodes and therefore do not optimize their transmit power 
values: p = Po for some pre-determined power values pg. If the S-D pairs aie non-cooperative, 
then each S transmits at full power and thus p = IP™"^. Given the transmit power values of 
the source nodes, we maximize the achievable rate by choosing the relax matrix. The PB problems 
are formulated into quadratically constrained quadratic programs (QCQP) and are then relaxed to 
semi-definite programs (SDP). The relaxed problems are convex and can be solved using efficient 
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convex optimization tools such as CVX [32]. We show that in some scenarios, the optimal solution 
of the relaxed problems attains the optimality of the original problems. In other words, the convex 
optimization methods solve the original problem efficiently in such scenarios (Please refer to Corollary 
[T] and Lemma |4] for details). The procedures used to obtain the PB are summarized in Algorithm [1] 

Algorithm 1 The pseudo-code for the Pareto boundary optimization in ([T6l) and (fTSl) . 
1: for j = 2 ^ K do 

\h ■ Rg PP- 

2: Compute the single-user-point of user j: j'^'^^ = maxR "|| g^ R [p'+i — ~' which the user 
j is the only user in the system with no interference. 



4: end for 

5: Define a tuple s to be a vector of possible values of 7^, s = [72, • • • , 7ii'] and s E § = 

V2 X V3 X • • • X Vi^. 
6: for each s G S do 

7: With input parameter s, matrix S is defined as in Solve the optimization problems (l22l) 
for general relay processing matrix optimization or (l24l for relay processing matrix optimization 
with interference neutralization. 

8: if the optimization problem is feasible then 

9: the optimal value 7I and s form a point on the Pareto boundary, which h a K — 1 

dimension hyper-surface, in the K dimensional space. 
10: else 

11: the values of s are unachievable subject to the constraints. 

12: end if 
13: end for 



A. The Pareto boundary: general relay optimization 

By introducing an auxiliary variable, we formulate the PB optimization problem with general relay 
processing matrix as a QCQP The optimization variable is of dimension + 1 as compared to 
number of elements in the relay matrix. Nevertheless, this provides a more structural formulation and 
amends the analysis as shown in the sequel. 

Lemma 2: The Pareto boundary of an IRC with instantaneous AF relay ([T6] ) is a rate tuple 

jj jj — H — 

(C'(7i ), (^(72), . . . , C(7i^)) in which 7{^ = Y^-^" Y and v is the optimal solution of the following 
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For a predefined integer N, let Vj = < 0, 
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optimization problem. The values 7^ , j 
problem. 



2, . . . ,K contribute to the constraints of the optimization 



max 



s.t. 



tH 



X12 V 



> 



7j, J 



(21) 



The matrices, for i = 1, 



X3 V < 0. 
, K, are given by 



X 



il 



(gri^gir)*(gri^g*r 



(gr 



X 



42 



K 

E 



^®g*rfhi. 
I "-jj I 

(gW®g?r)'' hil 



^M®{gir gfr) 
0ixM2 



1 



x. 



Q^^I 







lxA/2 



M2 X 1 



-Pr 



Proof: See Appendix IVII 
Define V = v v^, we obtain the following convex problem after removing the constraint rank(V) 



1: 



max tr (Xn V) 



(22) 



s.t. tr(Xi2V) = l, 

tr((X,i-7,X,-2)V)>0, j = 2, 
tr (X3 V) < 

which is a semi-definite program (SDP) as matrices Xii,Xj2,X3 are Hermitian and can be solved 
efficiently using SDP solvers. 

Corollary 1: By [33, Theorem 3.2] the rank of the optimum solution of (l22l ) is smaller than 
+ I. In the scenario of K = 2 S-D pairs, the rank of the optimal solution is one which means 
that the relaxation is tight and the obtained solution is the global optimal solution of (|2T]) . In occasions 
when the optimal solution returned by CVX is not rank-one, one can find a vector which satisfies the 
same constraint and objective values in (|22] | using the rank-one reduction procedures [34, Theorem 
2.3]. Such vector is thus one of the global optimal solutions of (|2TI) . For more S-D pairs, the result of 
SDP in (I22I ) is not rank-one. However, one can apply randomization approximation techniques [35] 
to approximate the optimal solution of (1211) . 
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B. The Pareto boundary with IN 

We manipulate the problem ([TST i in the same fashion as in (|2TI) . The optimization problem ([TST i has 
one more constraint (the IN constraint) and is thus optimized in a smaller dimension + 1 rather 
than + 1. Note that the condition given in Theorem [T] is a necessary and sufficient condition for 
the feasibility of IN (non-empty feasible set in (ITSl)). In other words, if the condition is not satisfied, 
then the optimization problem in (|23] | is not feasible regardless of the target SINR values 72, ... , ^k- 

Lemma 3: The Pareto boundary of an IRC with instantaneous AF relay and IN (fTSl) is a rate tuple 
(C(7*), (7(72), . . . , C{jk)) in which 7* = y^Bjiy and y is the optimal solution of the following 
optimization problem. The values 7j,j = 2,. . . ,K contribute to the constraints of the optimization 



problem. 



max 

s.t. 



y^Bny 



y^Bi2y = 1, 
y^B. y >0, 

y^Dsy <o, 

y^D4y = 0. 



j = 2,...,K, 



(23) 



where 



Bii 



L"^eie^|"L L-^eihn 



till 



Pi, 



B 



12 



1 



B. 



e, ej L -7, ( G,,* GJ » ej 



T 



7j 



Q 0i^2xi 

Opmax 

T^T T^Tvec(H) 
vec(H)^ T vec(H) T vec(H) 
Proof: See Appendix IVIII-CI ■ 
Let Y = yy^. Using SDR techniques, the problem in (1231 ) can be relaxed to the following 
problem by dropping the rank one constraint on Y. Problem (l24l) is a convex problem, in particular 
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a semi-definite program, which can be solved efficiently. 

tr fBiiY' 



max 

s.t. tr(Bi2Y)=l, 

tr(B,Y)>0, j = 2,...,K, (24) 
tr^DgY) <0, 
tr (p^ y) = 0. 

Note that, if the optimal solution in (l24b is rank one, then such solution solves (1231 ) optimally. If the 
optimal solution in (l24l) is not rank one, then the optimality of (1231 ) can no longer be guaranteed. In 
the following, we characterize the rank of the optimal solution of (l24l ). 

Lemma 4: The optimal solution of i^-user IRC Pareto boundary problem with IN in (|24] |. Y, 
satisfies 

rank(Y) < y/K+l. (25) 
Proof: Note that the matrix D4 ^ in (l24l) with rank — K and therefore we have eigenvalue 
decomposition of D4 as 

64 = [V, Vo] diag(Ai, . . . , Xk-~k, ^ix({k-+i)-(k-~k))) [V^; V^] (26) 

where V is the eigenvector matrix of D4 corresponding to the eigenvalues Ai > . . . > Xk'^-k > 0; Vq 
therefore spans the null space of D4. Let y = Vq x for x G ([^{K'^+'^)x'i- ^j^j rewrite the optimization 
problem to 



max 



x^V^BnVox 



s.t. x^Vf Bi2Vox = 1, 



(27) 



x^ (b,i - 7,B,-2) Vo X > 0, j = 2, . . . , 
x^VoD3Vox<0. 

There are now K + 1 constraints in (|27] ) and result follows from [33, Theorem 3.2]. ■ 
Corollary 2: The optimization methods and results for (|24] | are similar to Corollary [T] and shall 
not be repeated here. 



V. The optimal transmit power at informed simple S-D nodes 

In the scenaiio where the source and destination nodes are infomied of the presence of the relay 
and are willing to improve the rate performance of the system by cooperation, we can improve further 
the Pareto boundary by optimizing the transmit power at the source nodes. In the i^-user SISO-IC, 
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the power allocation at the PB is characterized and is obtained by searching over K — 1 real-valued 
parameters [31]. In the following, we observe that given the relay matrix R, the Pareto optimal source 
power in IC does not apply to the IRC. The is due to the dependence of the source power in the 
relay power constraint. Nevertheless, we obtain the Pareto optimal source power as a function of 
relay matrix R and target SINR values ^2, - ■ ■ ,1k- 

Theorem 2: For any given relay matrix R, the PB of (fT6l ) is a set of rate tuple (C(7*), C(72), . . . , C{ik)) 
which are attained by the optimal transmit power at S nodes, p# = [pf ,'P2-k\ — ^' 

' .* min (r-- ~ ''^'^'^''^ ~ ^ min ( ""^^^ ^^-'^""^'^ " ^^^^ ^ 

r ' (ll..P-c^[A]^,W)a,) '"^^^n [^]- 

p*^ = min {[A]^Xk) (q - Pf) ' Pr 



(28) 



where = [|| gr.2 f , • • ■ , II grX f ] - = 7m+i (llg(L+i)rI^II^ + 1) and 



[AU={ .W)(n.+i)+g&+i).Rg.(™-M)l^ ifm + 1-/, ^29) 

-7m+l|/i(m+l)i + g(^+i)r Rgri P if (m + l)^ I. 

and [A](2:ft:) is a matrix which consists of the second to K-th columns of the matrix A. 

Proof: See Appendix IVIII-Dl ■ 
Theorem |2] gives the Pareto optimal transmit power p* as a function of the given relay matrix R and 
the target SINR values 72, ... , ik- Similar results can be obtained when IN is employed, as shown 
below. 

Theorem 3: For any given relay matrix R, the PB of (fTSl ) is a set of rate tuple (C(7*), (7(72), . . . , CdK)) 
which are attained by the optimal transmit power at S nodes, u = [ui, . . . , uk] > 0: 



pmax _ ^i^(s) R(s)^) ^ II g^j R(s^ll2 



' l|gr.lR(s)P ^ l|grlR(s^"' ' ' 

7,(||gf,R(s)f + 1 



max 



(30) 



Proof: See Appendix IVIII-El 
Note that Theorems |2] and |3] give the Pareto optimal source power allocation for a given relay 
processing matrix R. While the joint optimization of the source power and relay processing matrix 
is highly complicated, one can approach the problem by solving iteratively (a) the relay processing 
matrix given the source power (l22i and (l24l and (b) the source power given the relay processing 
matrix (|28] ) and (|30] |. However, the iterative optimization approach may not converge to the global 
optimal solution. 
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In the following section, we provide numerical evidence of performance gain of a relay introduction 
to a SISO interference channel. In particular, in a setting of uninformed source nodes, we show the 
rate improvement of solely introducing and optimizing the relay strategy whereas in a setting of 
informed S-D nodes, we compare the rate performance of the general relay optimization and the 
relay optimization with IN to the rate region of a SISO IC. 

VI. Simulation Results 

For illustrative purposes, we let K = M = 2. We assume that each element in the channel matrices 
H, G-rt, Gdr is an independent identically distributed complex Gaussian variable with zero mean and 
unit variance. In Section IVI-AI we simulate the achievable rates of the SISO IC (marked as squares). 
For comparison, we simulate the achievable rates of the fully connected IRC with the same S-D nodes, 
by introducing a relay, equipped with 2 antennas, to the aforementioned IC, and the relay can choose 
to enforce IN (marked with asterisks) or not (marked with triangles). We show that optimized relay 
strategies improve achievable rate regions. In Section IVI-BI and IVI-CI we compare the average sum 
rate and proportional fairness utility achieved by optimized relay strategies and by power allocation 
on the IC respectively. In Section IVI-D[ we illustrate the sum rate performance and proportional 
fairness utility when the transmit power constraint at source nodes and relay nodes vary. 

A. Rate region improvement 

In Fig. |2l we plot the achievable rate region of a two-user SISO IC with transmit power constraint 
at each source node P"^"-^ = IQdB. Introducing an instantaneous relay, equipped with 2 antennas, we 
obtain an IRC. We set the relay power constraint as P™"^ = 20dB. The achievable rates achieved by 
general relay optimization and IN outperform the IC case. The black arrows originate from the Nash 
Equilibrium point: the rate points in which both users transmit with full power. The north-east side of 
the arrows mark the rate region improved by the relay, in the scenario of uninformed source nodes. 
This validates our intuition that optimized relay strategies can improve achievable rates of the system 
even if the source nodes are oblivious to the existence of the relay and do not change their transmit 
power. Further note that the single user points achieved in IRC with general relay optimization always 
outperform the single user points in a SISO IC. It demonstrates that the relay not only is capable of 
reducing interference in the system but also forwarding the desired signal to the destinations. 

In Fig. [21 we reduce the relay transmit power to P^"-^ = WdB. We observe that the rate region 
achieved by IN reduce significantly because the relay is not able to neutralize interference and improve 
desired signal quality with limited power. 
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Fig. 2. The rate improvement of relay optimization on a two-user SISO IRC with K = M ^ 2, P^""^ = 2MB, 
pmax _ iQfj^]^ jYie arrow marks the increment of rate region by introducing a relay into the system and optimizing the 
relay strategy. 




Fig. 3. The rate improvement of relay optimization on a two-user SISO IRC with K = M ^ 2, P™""" = lOdB, 
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B. Average sum rate improvement 

In Fig. m we show the maximum sum rate achieved by general relay optimization, IN and power 
allocation on the IC, averaged over 100 independent channel realizations. The power constraint at 
the source node is assumed to be p™'*^ = \{)dB and we increase the relay transmit power from 
bdB to 25dB. We observe that the optimized relay strategy without IN always outperform the 
maximum sum rate of the IC, demonstrating that an instantaneous smart relay can improve average 
sum rate performance. Further, we observe that the performance of IN is limited by the relay transmit 
power. Although IN is analytically appealing, there are limitations of the implementation of IN. Such 
scenarios include strong interference channels in which the receivers have strong interference from 
other transmitters in the system. In this case, more power at the relay may be required to completely 
null out interference and if such power is not available to the relay, then IN is not feasible. On the 
other hand, if the strength of the interference channel is not strong, enforcing IN, the relay loses its 
optimization degrees of freedom and may not be able to achieve some operating points as the general 
relay optimization would achieve. 




Fig. 4. The average sum rate a two-user SISO IRC with K = M ^ 2, pj""^ = IQdB. The optimized relay strategies 
improve the average sum rate of the system significantly. 
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C. Proportional fairness improvement 

While average sum rate is an important system performance measure, user fairness holds importance 
in many applications. In Fig. [51 we illustrate the average proportional fairness utility which is defined 
as the maxijj /{^(i?! — R^^){R2 — RJ^""^). We observe that the optimized relay strategies, with and 
without IN, provide promising proportional fairness and better sum rate performance compared to 
IC. 




D. Performance measures in terms of transmit power constraints 

It is interesting to observe that the performance of optimized relay strategies depend on both P"^"-^ 
and P^""^. This is due to the amplify-and-forward nature of the relay. If the transmit power from 
source nodes is high, the relay can spend less power on amplification of signals due to the relay 
power constraint. In Fig. |6l we plot the maximum sum rate achieved by optimized relay strategy with 
and without IN and the maximum sum rate in IC. Note that the feasibility conditions of IN, shown 
in Theorem [U is validated in Fig. |6] For a fixed relay power i-"™"^, when the source power increases 
such that the conditions are violated, IN is not feasible and the achievable rate is zero. For the general 
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relay optimization, for a fixed relay power and increasing source power, the rate performance is not 
always increasing because the increased source power increases the interference power and the relay 
may not have enough power to manage interference and amplify desired signals simultaneously. 

In Fig.|7] we show the maximum proportional fairness utility achieved by optimized relay strategies 
and IC. When both source power and relay power are abundant, the fairness is desirable. However, 
when the source power (which is also the strength of interference) is relatively stronger than the relay 
power, the fairness achieved by the relay strategies is overtaken by the proportional fairness utility 
achieved by IC. 



H IRC with IN 
niRCw/o IN 
IC 




Fig. 6. The sumrate of a particular channel realization of a two-user SISO IRC with K = M = 2. 



VII. Conclusion and Future research directions 

The achievable rate region of a SISO-IC has been an on-going research topic, with recent interest 
on the question whether a relay introduction to the SISO-IC, obtaining an interference relay channel, 
provides any performance gain. In this paper, we study this problem by assuming an instantaneous 
amplify-and-forward relay with uninformed source and destination nodes in the system. We examine 
the gain of rate region of the relay introduction by formulating the Pareto boundary problem with 
optimization over relay processing matrix. The optimization problems, with and without the employ- 
ment of interference neutralization techniques, are solved using semi-definite relaxation techniques. 
The global optimality of the solutions are proved in the scenario of two source and two destination 
nodes. In the scenario of informed source nodes, we allow the source nodes to optimize their transmit 
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Fig. 7. The proportional fairness utility {Ri — R^^){R2 — R2^) of a particular cliannel realization of a two-user SISO 
IRC with K = M = 2. 



power. The transmit power values at the source nodes which attain the Pareto boundary are obtained 
in closed-form. Simulation results confirm that instantaneous relay is able to improve the achievable 
rate region, even in the scenario of uninformed source nodes ; improve average sum rate and average 
proportional fairness of the system. 

This paper motivates the study of performance of the IRC with an AF relay which can be 
implemented easily in practical applications as the relay is only responsible for a simple forward 
process which does not incur a processing delay as compared to other complicated relays. As a 
preliminary study, we only allow the relay to choose between IN or no IN. To evaluate the full 
potential of the AF relay, one may allow the relay to choose different relay strategies, e.g. interference 
forwarding, signal amplifying, interference neutralization, etc., depending on its power budget and 
the channel qualities in the system. Another interesting problem is the physical placement of the relay 
with the goal of rate performance improvement. 

VIII. Appendix 

A. Proof of Theoren{J\ 

The goal of this section is to prove that the necessary and sufficient condition of S, satisfying 

vec(S)^Qvec(S) < f,"^"^ 

< 

Tvec(S) = -Tvec(H) 
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with Q = (diag(p) + G-,^) (g,; G^,^), is 

(Tvec(H))^ (TQ"^T^)"^Tvec(H) < (31) 

Perform eigenvalue-value decomposition on the Hermitian matrix Q = T U and we let x = 

r^/^ u vec(S) and x < fT*"^'. Let F = r^^/^ ^nd we have 

Tvec(S) = TFx = -Tvec(H). (32) 

Note that the matrix T F is of dimension K{K — 1) x K"^ and has nullity of K. Denote the null 
space of T F by J\f{T F) and we choose 

x = x„ + x/„ x„ G AA(TF), x;, ^ AA(TF). (33) 

From (l32l ) and (|33] ). we have 

x;, = -(TF)tTvec(H). (34) 

The sufficient condition of IN is thus 

||x,,||2 = ||(TF)tTvec(H)f = (Tvec(H))^ (TF(TF)^)"^ Tvec(H) < P,'^'^^ (35) 

Once x/i is chosen according to (l34b and satisfy (|35] |. we are free to choose Xn G 7V(T F) as long 
as II xf < P^^"^. 

To see that (l35l) is a necessary condition, we need to prove that if IN is feasible then (l35l) must be 
true. We prove by contradiction and assume that IN is feasible and it is possible to find a solution x 
such that II X f < P™"^^ but (Tvec(H))^ (TF(TF)^)"^ Tvec(H) > P™"^. By assumption, there 
exist X such that -Tvec(H) = TFx. We multiply both sides by (Tvec(H))^ (TF(TF)'^)"^ 
and we have 

(Tvec(H))^ (TF(TF)-f^)"^Tvec(H) -(Tvec(H))^ (TF(TF)-f^)"^ TFx > P™"^. (36) 

Since the product on the left of the equality (a) is a real number, the product on the right side of 
(a) must be a real number also. Thus, we can write 

x = -a(TF)-^ (TF(TF)^)"^Tvec(H), ,aGM+. (37) 

Substitute (|37] ) into (l36l ). we have 

a(Tvec(H))^ (TF(TF)-f^)"^ Tvec(H) > P™'^^'. (38) 

By assumption, (Tvec(H))^ (TF(TF)'f^)~^ Tvec(H) > P;""^ and thus a > 1. Computing the 
norm of x, we have ||xf = a2(Tvec(H))^(TF(TF)^)~i Tvec(H) > P;"'^^ which violates 
the power constraint. Thus, (Tvec(H))'f^ (TF(TF)^)"^ Tvec(H) < P™"^ is a necessary and 
sufficient condition for interference neutralization. 
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B. Proof of Lemma |2] 

In this section, we show that steps to obtain (|2T]) from ([T]). From (|3]), the desired signal power of 



(39) 



which is due to a B c = (c(g) a) vec(B). The interference and noise power of Tx-Rx pair i is then 

K 

Y.\hii + (Sh ® ^iT vec(R)|2p; + II g^ R ||2 + 1 

(40) 

K ^ ' 

= J2\hii + {gri » Stf vec(R)|2p, + II vec(R)^ (g*, ^ Ik) f + 1 

where the manipulation (a) is due to g^R = vec(R^)^(g*^ (g) 1/^) = vec(R)^ Z^(g*^ (g) 1/^) and 
the matrix Z which satisfies vec(R^) = Zvec(R) is called a commutation matrix [28, Section 9.2]. 
Z has elements zeros and ones and Z = Z^. Utilizing the Chames-Cooper Transform [36]-[38], we 
let r]/t = vec(R) where t € C and v = [r]'^,tf e c(^'+^)'<^. The SINR of user i is 

\hiit + {gri'^^S*rfv\^Pi 



SINR, 



E^i \hat + (gw €D g*,)^ V \'Pi + II (g; » Ix) P + |tP 

(^) \hiit + {gri^g*rfv\'^Pi 

Eg. |/i.;i + (gw ^5 g*,)^ r/ |2P; + r?^^ (I^ g,, g^) r? +|t|2 

X,;i V 



(41) 



X,;2V 

where (a) is due to Z ( A (g) B) Z = B A and 



X 



il 



(gri g*r)* (g« g*irf (g« ® glr)^ ^« 



K 



iSrl gir)* (gr/ S*irf iSrl ^ S*irf hu 



Pl + 



lM^(girg*>) Oa/ 



xl 



0lxM2 



(42) 



The power constraint at the relay is 

tr (R Q R^) = vec(R)^ (Q^ 1) vec(R) < P; 
(Q^(8)Im) Ojvpxi 

From Eqt. (ITO and (l43T l. we obtain the formulation in Lemma |2] 



(43) 



V < 
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C. Proof of Lemma\3\ 

In this section, we show the steps to obtain (1231 ) from ([TSl l. We rewrite the noise power to the 
following 

R = g^ {g'^^ s g^/^ (g;^;^ s g^/^ g^^ 
- „H (^-H ar-^ r-^ r-^ a- 

vec(S^)^ (G^; g*, G,-/) (G,-; g*, G,"/)^ vec(S^)* 

veciSf Z (g,; g;, gl G-J G,-/ G;,^) Z vec(S)* (44) 



vec 



(S)^ (g;/ G;,^ G^; gl gl G^J) vec(S) 



vec(S)^ (g,-; G;,^ ® G,; g„ gl G,,^) vec(S) 

= vec(S)^ (g-* G-J ® ei ef ) vec(S). 

The transition (a) is due to the Kronecker product properties, a^B^C = vec(B)^ (a C) [39]. The 
commutation matrix Z of dimension K x K satisfies Zvec(S) = vec(S^), [28, Section 9.2]. The 
transition (6) uses such properties and (Ai (g) Bi)(A2 (8> B2) = Ai A2(X'BiB2. Then in transition 
(c), we use the property Z(A(g)B)Z = B (g) A. In transition (d), we uses the fact that the noise 
energy is a real scalar and a complex conjugation does not affect its value. The last equality is due 
to that fact that gj^ is the i-th column of G^r and G'^^ g^^ = [G^^^ G^^] (-. ^-^ = e^. 

Now, we rewrite the SINR and power constraints as a function of vec(S). The signal power at Dj 
is rewritten as 

\hii + Si\'^ = + ef Lvec(S)|2Pi (45) 
where L is a row selection matrix, s = Lvec(S). From (|44] | and (l45T l. the SINR of Di is 

vec(b)^ Bj2vec(S) + 1 
Recall from (llScI) and (llSdl ) that any feasible solution satisfies 

vec(S(s))^Qvec(S(s)) < P,"^'^^ 

(47) 

T vec(S) = -Tvec(H) = b. 
D4 

However, the second constraint (1471 ) creates complication to the optimization problem because of 
the asymmetric structure of D4. We here propose an equivalent constraint 

||D4vec(S) -bf = 0. (48) 
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From Eqt. (|46] | and Wfl and let w = vec(S), (fTSl ) is equivalent to the following formulation, 

|/iii + ef Lw pPi 



max 

S.t. 



^Bi2W+l 



\hjj + ejLw\^Pj 



w 



w < PI' 



^B.2W+1 



> 



7j, 



j = 2,...,K, 



(49) 



II D4W-b II = 0. 

In the following, we convert the optimization problem in (l49l ) into the standard QCQP. We proceed 
with the Chames-Cooper transform [36]-[38] by substituting the optimization variable w = wt for 
some complex scalar t / 0. The optimization problem is rewritten as 

[efL,/in]^[efL,/in]Piw 



rnax 



s.t. 



W^Bi2 W + t2 



w 



H 





H 






ej L, hjj 




ejL, hjj 


PjW 



Bj2 W + t2 



> 



7j> 



(50) 



II D4W - btf = 0. 

We denote Bji = [e^ L, /ijj] ^ [e?^ L, /ijj] Pj, i = l,...,K. Without loss of generality, set the 
denominator of the objective function to one and define a new optimization parameter y = [w^,t]-^. 
We obtain (l23]l with 



B 



il 



h*efL 



\h. 



"I '■'■It 

|2 



B 



i2 



B, 



Bi2 Oa'2 X 1 

Olxii-2 1 



T ^ T 



Olxifa 1 

T ^.^Tt _^,.( n-* n-T ^ „ . „T 



D3 

Oixi^^ 
Df D4 



0_ff2xl 



e,- ej L -7,- (^G,/ G,/ e,- ej 



ej hjj- 



(51) 



Df b 



b^B4 



b^b 



Tvec(H) 

vec(H)^ T vec(H) T vec(H) 
Note that all above matrices above are Hermitian matrices. 



0lxK2 



Pt + 1 



25 



D. Proof of Theorem |2] 

For any given relay matrix R, we write the optimization in ([T6] l as 

pei?- E/li,^^i + gf. Rgw + II gf. R P + 1 



max ^ — „ ^„„, (52a) 



s.t. |/^,,+gf,Rg.,pP, ^ j = 2,...,i^, (52b) 

E^M^,>.7 + gf.RgwPP. + l|gf.Rf + i" ' 

tr I^R [jl^ g,, g^[ Pi + R^^ < Pr", (52c) 

p^<praax^ / = l,...,i^. (52d) 

If the problem is feasible, at optimality all the SINR constraints are settled in equality. To see this, 
denote the optimal power allocation as p and assume that SINR2(p) > 72 and SINRj(p) = 7j 
for j = 3, . . . jiT. Note that SINR2(p) is monotonically increasing with P2 whereas SINRi(p) and 
SINRj (p) are monotonically decreasing with P2- Thus, the decreased value of P2 increases both 
SINRi and SINRj , j = ?,,..., K. On the other hand, the power constraints (I52ct and (I52dt and the 
SINR2 constraint remain valid. This contradicts to the assumption that P2 attains the optimal point. 

Since all SINR constraints are active at optimality, we write all — 1 constraints (I52bl ) in the 
following: 

Ap = q (53) 

where for m = 1, . . . , K — 1, I = 1, . . . ,K 

^(m+l)(m+l) + g(L+l)r ^^r{m+l) P if m + 1 = /, 

-7m+l|^(m+l)i + g(L+l)r ^^rl ? if {m + I) ^ I. (54) 



Hm = 7m+l (ll g(L+l)r ^ f + l) • 

Note that the matrix A has dimension {K — 1) x K and we denote the i-th column of A as and 
the second to last elements of p as P2:A'- We have 



P2:K = [a2, . . . , ax] [q,- ai A j • (55) 

For the brevity of notations, we let [A]^. ^-k) — [^2t ■ ■ j^k]- Substitute into the power constraint 
(I52cl ) and denote c-^ = [|| g^2 IP) • • • > II gr/r IP] ; we have 



g,i R f A + [A]^;2^^) (q - ai A) < Pr'' - tr (RR^) 

prnax _ tr (R R^) - [A] ^^^2:/^) Q (^6) 



^ Pi < 



g^RP -C^ [A](^_2:/0^1 
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From (|55] ) and (I52dl i. we have 

P2:K = [A]^U) (q-ai A) < Pr"'l(i^-l)xl 

^ ^ (57) 

4^aiA>q-Pr" [A](,2:i^)l(i,-i)xi 

Note that ai = [-72I/121 + Rg^i P, ■ ■ ■ , -7i<-|/ii^i + gf^RgriP]^ and we have K-1 upper 
bounds of Pi and for 1 < /c < — 1: 

Pi < - — ^^rf ^ ' . (58) 

On the other hand, the objective function SINRi is monotonically increasing with Pi- To see this 
mathematically, we write 

l/tii+ggRg^iPA 

\hii + gf. Rgw PA + II gf, R P + 1 

^ l^ii + ggRg^iPA 

b^P2:i^ + l|gf.RP + l 

|/.ii+gf,Rg,i|2A (59) 



(q-aiA) + ||gf,RP + l 
l/iii+gf^Rg^iPA 



- A + (b^ [A]iXk)^ + \\ gf. R P + 1^ 

where = [|/ii2 + gf^ Rgr2 • • • , \hiK + g(^. Rg^i^ ■ The last equality is of the form ^l^^ 
which can be manipulated as 

Zl A Zi Z2P1 + Zs- Z3 Zi f ^ _ ^3 \ ^^Q^ 



2;2Pl + Z3 ^2 Z2P1 + Z3 2:2 V Z2P1 + Z^, 

which is indeed monotonically increasing with A- From (l56l ) and (l58T l. we have: 

g,iRP-c^[A]^J^^)ai^ ^=i.-'^V [A].i 



(61) 



E. Proof of Theorem \3\ 

With the requirement of IN, the interference from Tx j at Rx i is canceled and therefore the 
function SINRj(s, Pj) is independent to the transmit power from any other transmitters j ^ i. From 
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([T8] ). for any given s, we have 

|/iii +gf,R(s)g^ip P] 
gf(.R(s)P + l 



+ 



+gf^R(s)g„-|2p- 
s.t. --„^.- - >7. .-2,...,^, (62b) 

tr (n g^ Pi + R^^ < P--, (62c) 

p^^pniax^ / = !,... if. (62d) 

The IN constraint disappears from the above optimization problem because it is independent to p. 
Denote the optimal power allocation by u = [ui, . . . ,uk\- We can observe that SINR constraints 
(I62bl ) and power constraint (I62cl i must be active at optimality. Otherwise, let SINR2 > 72- The value 
of U2 can be decreased by a very small amount e without violating (I62bl ) and ui can be increased by 
e without violating (I62cl ). This new increases the objective value and leads to contradiction that 
we are at the optimal point. From (I62bl i and (I62cl i. we have 



7,(||gf,R(s)f + 1 



(63) 



5] II g,, R ||\i = - tr (R(s) R(s)^) . 



Therefore, 



. I Pr"-tr(R(s)R(s)^) ^ ||g,,R(s)f , 
Ml = mm L^, — > 71 — . ,,„ Ui,P' . (64) 

I l|grlR(s)P ^^||g,iR(s)P ^' « ; 
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